Neutron to Mirror-Neutron Oscillations in the Presence of Mirror Magnetic Fields 



I. Altarev,! C. A. Baker,^ G. Ban,^ K. Bodek,-* M. Daum,5'El P. Fierlinger,*^ P. Geltenbort,"^ 
K. Green,^'^ M. G. D. van der Grinten,^'* E. Gutsmiedl^ P. G. Harris,* R. Henneck^ M. Horras,'''^ 
P. Iaydjiev,2'[t| S. Ivanov,^'! N. Khomutov ^ K. Kirch,5.| S. Kistryn,^ A. Knecht,^'^ P. Knowles/° 
A. Kozela," F. Kuchler,^ M. Kuzniak,^' ^'B T. Lauer/^ B. Lauss,^ T. Lefort,^ A. Mtchedlishvili,^ 
O. Naviliat-Cuncic,^ S. Paul,^ A. Pazgalev,^° J. M. Pendlebury,* G. Petzoldt,^ E. Pierre,^' ^ 

C. Plonka-Spehr,i2 G. Quemener,!^ D. Rebreyend,^^ S. Roccia/^ G. Rogel/'^ N. Severijns," 

D. Shiers,* Yu. Sobolev,^^ R. Steepler, ^ A. Weis,^" J. Zejma,'' J. Zeniier,^^ and G. Zsigmond^ 

^ Technische Universitdt Munchen, D-85748 Garching, Germany 
'^Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX, United Kingdom 
''LPC Caen, ENSICAEN, Umversite de Caen, CNRS/IN2P3, F~14050 Caen, France 
''Marian Smoluchowski Institute of Physics, Jagiellonian University, 30-059 Cracow, Poland 
•''Paul Scherrer Institut (PSI), CH-5232 Villigen PSI, Switzerland 
^Excellence Cluster 'Universe', Technische Universitdt Miinchen, D-85748 Garching, Germany 
''institut Laue-Langevin, F-38042 Grenoble Cedex, France 
^Department of Physics and Astronomy, University of Sussex, Falmer, Brighton BNl 9QH, United Kingdom 

^ JINR, 14 1980 Dubna, Moscow region, Russia 
'"University of Fribourg, CH-llOO, Fribourg, Switzerland 
Henryk Niedwodmczanski Institute for Nuclear Physics, 31-342 Cracow, Poland 
'^Institut fiir Kernchemie, Johannes-Gutenberg-Universitdt, D-55128 Mainz, Germany 
'■^LPSG, Umversite Joseph Fourier Grenoble 1, CNRS/IN2P3, Institut 
National Polytechnique de Grenoble 53, F-38026 Grenoble Cedex, France 
'"^ Katholieke Universiteit Leuven, Instituut voor Kern- en Stralingsfysica, Celestijnenlaan 200D, B-3001 Leuven, Belgium 
'^Institut fiir Physik, Johannes-Gutenberg-Universitdt, D-55128 Mainz, Germany 

(Dated: August 18, 2009) 

We performed ultracold neutron (UCN) storage measurements to search for additional losses due 
to neutron (n) to mirror- neutron (n') oscillations as a function of an applied magnetic field B. In 
the presence of a mirror magnetic field B' , UCN losses would be maximal for B B' . We did 
not observe any indication for nn' oscillations and placed a lower limit on the oscillation time of 
r„„/ > 12.0 s at 95% C.L. for any B' between and 12.5 ^iT. 

PACS numbers: 14.80.-j, ll.30.Er, 11.30.Fs, 14.20.Dh 



I. INTRODUCTION 

The idea of restoring global parity symmetry by intro- 
ducing mirror particles dates back to Lee and Yang ,11]. 
In 0] , this idea has been significantly expanded and was 
later adapted to the framework of the Standard Model of 
particle physics A recent review can be found in 
Interactions between ordinary and mirror particles are 
possible, e.g., they both feel gravity, making mirror mat- 
ter a viable candidate for dark matter d, H, 0, 9 ■ Besides 
gravity, new interactions could lead to mixings between 
neutral particles and their mirror partners. 

Fast nn' oscillations were introduced in Q to explain 
the existence of ultra-high energy cosmic rays, based on 
a crude limit on the oscillation time t„„' ^ 1 s. This 
weak limit was one of the motivations to perform a first 
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dedicated measurement, which resulted in a lower limit 
of T„„, > 103 s (95% C.L.) [3. The experiment relied on 
comparing the numbers of stored UCN remaining after 
a certain storage time for zero magnetic field and for an 
applied magnetic field B of several /iT [11]. Only for B sa 
would the ordinary and mirror state be degenerate and 
nn' oscillations could occur leading to an additional loss 
of stored UCN. Shortly thereafter, an improved result of 
Tnn' > 414 s (90% C.L.) was report ed [13] and further 
improved to t„„- > 448 s (90% C.L.) p|~^ 

So far, the limits were obtained assuming a negligible 
mirror magnetic field B' , except from an attempt in [T3| 
for mirror magnetic fields in the range to 1.2 /iT. Here, 
we report the first systematic search for nn' oscillations 
allowing for the presence of B' . The basic measurement 
principle remains unchanged with the exception of scan- 
ning B in order to find a resonance of maximal UCN 
losses dX B K, B' instead oi B « 0. The limits on B' 
from, e.g., a limit on the amount of mirror matter in- 
side the earth [T^ are very weak. Photon-mirror-photon 
mixings could possibly provide an efficient mechanism 
to capture mirror matter in the earth allowing for B' of 
several [iB] ■ Mirror magnetic fields not bound to the 
earth are also conceivable and would additionally lead to 
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daily modulations in the UCN counts - an unmistakable 
signature of a possible origin of B' . In the following, we 
will first introduce the theory of nn' oscillations in the 
presence of B', describe the measurements and conclude 
with the two analyses conducted: i) the search for daily 
modulations and ii) the search for a resonance. 



II. nn' OSCILLATIONS IN THE PRESENCE OF 
A MIRROR MAGNETIC FIELD 

For the calculation of the nn! oscillation probability 
with finite B' , we follow the arguments of [l^. Defining 
2hu3 = finB and 2huj' = /i„-B' and introducing the oscil- 
lation time Tnn' and the Pauli matrices cr, the transition 
from the ordinary to the mirror state (and vice versa) is 
described by the interaction hamiltonian 



Defining a coordinate system with b = (0, 0, 6), b 
uj'\, and a = {ax,0, Oz), ax = 2\uj x + a;' 

(w^ — a;'^)/|a; + a;'|, leads to the 4x4 matrix 
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Ti can be diagonalised using a transformation ma- 
trix with mixing angles fulfilling tan 20 = l/(azT„„'), 
tan 20 = ax/{b — a^), and tan2(/)' — ax/{b + dz) with 
ciz — az \/l + l/(a2T„„/)2 [l^ . The eigenvalues of H are 
±2a) and ±2uj' given by 2uj = ax sin 2(j) + (b — dz) cos 20 
and 2a)' = ax sin 20' + (b + dz) cos 20'. The time depen- 
dent probability for the transition from n to n' is then 
given by 



Pr, 



,{t) = sin2(26l) [cos^(0-0')sin^(t/r_) 

+ Sin2(0-0')sin2(t/r+)] , (3) 



where t± = \u)±u)'\~^ are the effective oscillation times. 
The oscillation probability depends on the magnitude of 
B and B' , the direction of B' given by the angle (3 relative 
to the up-direction of B (see below) , the oscillation time 
Tnn', and the time t. 

During the storage of UCN inside a chamber, the rele- 
vant time t is the free flight time t f between wall collisions 
in which the wave function is projected onto its pure n 
or n' state. The loss rate of UCN due to nn' oscillations 
is thus given as 
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(4) 



where denotes the collision frequency and (. . .)t, the 
averaging over the distribution of free flight times t f dur- 
ing the storage time tg. 



There are two distinct regions for the evaluation of the 
nn' oscillation probability. The first is the off-resonance 
region. From evaluations of Eq. this holds for \B — 
B'\ > 0.4 /iT. In this region, the time dependent terms 
in Eq. ([3]) oscillate quickly and average to 1/2 over the 
t f distribution. The loss rate is then expressed explicitly 
as 



1 B'^ + B'^ + 2B' B cos /3 2n^ 
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On-resonance, \B-B'\ < 0.4 ^T, the first term in Eq. ^ 
dominates for most of the parameter space. For that part 
of the parameter space, we have « 0' and, since t/T_ is 
small, sin^ {t/T~) ~ (t/T„)2. Therefore, we can replace t 

in Eq. ^ by ^{t'j-}t^ and write the loss rate as 
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The validity of Eq. ([6]) was checked by comparing to a 
full averaging over a realistic tf distribution. Deviations 
were less than 1%. Anyhow, our final limit is based on 
calculations using Eq. 

In order to obtain the values for (t/)t^ and \J {t^f)ts , 

a detailed Monte Carlo simulation of the experiment 
was performed using GEANT4UCN [16,] with parameters 
tuned to reproduce experimental data (such as character- 
istic time constants for filling, emptying, or storage). The 
t f distributions were obtained from the time of the refiec- 
tions of individual trajectories inside the storage cham- 
ber. Results are given in Table U for the two storage 
times ts used in the measurements. We varied the pa- 
rameters of the simulation in ranges still reproducing the 
experimental data to assess the systematic uncertainties. 
The number of surviving UCN after storage is 



N(tl)^N' exp{-RtX) 



(7) 



where Nq ^ is the initial number of UCN reduced by the 
usual losses during storage, and t* is the effective storage 
time for the UCN, including not only the time when the 
neutrons are fully confined, tg, but also the effects of 
storage chamber filling and emptying. The values for t* 
are given in Table IH 

In the case of a mirror magnetic field not bound to the 
earth, the observed neutron counts could be modulated 
with a period corresponding to a sidereal day {dsid — 
23.934 h) as the angle (3 would be modulated. For the 
off-resonance case, the observed counts are then given by 
N{t) = C + Aj^^ cos (27r(i - io)/4»d) with 
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TABLE I: Results for {tf)t, and yj{t])t, using Monte Carlo 

calculations and the effective storage times t^. The values 
at the right side of the arrow denote the values used in the 
calculations in order to obtain a conservative result. 



ts [s] 


75 


150 


{tj)t^ [s] 


0.0403(4) -» 


0.0407 0.0442(4) 0.0446 




0.0532(5) ^ 


0.0527 0.0586(6) 0.0580 




98(3) ^ 


95 173(3) ^ 170 



B' and B'^ are the components of B' parallel and per- 
pendicular to the earth's rotation axis, A the latitude at 
the experimental site, the phase, and the — (+) sign 
stands for magnetic field up (down). 



III. MEASUREMENTS 

The UCN storage experiments were conducted at the 
PF2-EDM beamline 'lg| at the Institut Laue-Langevin 
(ILL) using the apparatus for the search of the neu- 
tron electric dipole moment [l^. The main features of 
the apparatus are: i) the possibility to efficiently store 
UCN in vacuum in a chamber made from deuterated 
polystyrene [l9| and diamond-like carbon and ii) the sur- 
rounding 4-layer Mu-metal shield together with an inter- 
nal magnetic field coil that allowed to set and maintain 
magnetic fields with a precision of ^ 0.1 /xT. A typical 
measurement cycle consisted of filling unpolarised UCN 
for 40 s into the storage chamber of 21 litres, confining the 
UCN for 75 s (150 s) and subsequently counting ~ 38000 
(~ 24000) UCN over 40 s in a ^He detector tloj]. For a 
given magnetic field value, we always performed 8 cycles 
with a storage time of 75 s and 8 cycles with a storage 
time of 150 s. After these 16 cycles, the magnetic field 
direction was changed from up to down and measured 
again for 16 cycles. The averages of the different B field 
settings, applied randomly, were 0, 2.5 ^T, 5 /zT, 7.5 /iT, 
10 /iT, 12.5 /iT. Before doing a zero field measurement, 
the 4-layer magnetic shield was demagnetised resulting 
in i? < 50 nT. In total, data taken continuously over 
approximately 110 hours was used for the analysis. 



IV. NORMALISATION OF THE UCN DATA 

The data showed a trend to higher UCN counts over 
the course of the measurement period. The increase 
amounted to ~2.5% for 75 s storage time and ~5% for 
150 s storage. We attribute this increase to slowly im- 
proving vacuum conditions inside the chamber. A com- 
bined fit to both data sets was performed with the func- 
tion 

ftSt) = Nt^ exp [-Cptse-^ - Cflt^g-^) (9) 



TABLE II: Results of the fits using Eq. ([TOj to the up/down 
asymmetries A for the five different magnetic field values and 
the upper limits on the amplitude of a daily modulation 
at 95% C.L. and for any value of the phase to- 
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18.05/12 


5.0 


12.5 


1.0 ± 1.7 


17.1 ±9.8 


10.13/12 


5.0 



with two normalisation constants iV75 and iViso and two 
constants proportional to a decreasing overall pressure 
Cp (with a characteristic time Tp) and a decreasing out- 
gassing rate Cr (characteristic time tr) of the storage 
chamber, which is sealed off from the pumps during stor- 
age. The per degree of freedom, 1386/1204, is satis- 
factory. Assuming a UCN loss cross section per molecule 
of 0(10 b), the fitted constants Cp and Cr translate into 
an initial pressure of O(10~^ mbar) and an initial out- 
gassing rate of O(10^^mbarls~^ cm~^) which both seem 
realistic [l^ . We normalised the UCN counts for a given 
cycle by the prediction of Eq. ^ and slightly increased 
the statistical error by adding the fit error in quadra- 
ture. Residual drifts (< 0.5% over several hours) showed 
a weak correlation to the ILL reactor power. Their effect 
on the final result is negligible. 



V. ANALYSIS 

We conducted two different types of analyses: i) The 
search for a modulation in the UCN counts and ii) the 
search for a resonance in the UCN counts as a function 
of B. It is clear from Eqs. ([5]) and ([5]) that the res- 
onance analysis will always be sensitive to nn' oscilla- 
tions regardless of the origin of the mirror magnetic field 
and possible modulation periods whereas the modulation 
analysis is not. In Eq. ([5]), cos/3 will either be a fixed 
value or the average over a modulated cos [3. Addition- 
ally, the amplitude of the modulation tends to zero for 
small B' and the constant term C of the oscillation prob- 
ability is for all parameters larger or equal to the modu- 
lated part A {B'^ + B'^ > 2B'B cos 13). Given the same 
statistics and no systematic errors from averaging over 
longer periods, the resonance analysis will always yield 
tighter constraints on r„„' than the modulation analysis. 
As a means of crosschecking and discovering the possible 
origin of B' , both types of analyses have been performed. 



A. Search for a Daily Modulation 

In order to search for a modulation without be- 
ing affected by the slow residual drifts present in the 
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FIG. 1: (colour online) Combined fit to the normalised UCN 
counts as a function of applied magnetic field B for 75 s (dark 
green squares and solid line) and 150 s (light green triangles 
and dashed line). Positive (negative) B values correspond to 
B field up (down). 

normalised UCN data, we calculated the up/down- 
asymmetries in the UCN counts A = {N'f-Ni)/{N^-+Ni) 
from the two subsequent (within ~1 h) measurements at 
B field up and down. The two asymmetry data sets for 
75 s and 150 s were separately normalised in order to 
have zero weighted means. A modulation in the UCN 
counts would show up in the asymmetry with the same 
amplitude A as given in Eq. ((HI): 



A{t) = ^77^ cos [^{t~ to) 



(10) 



We searched for a modulation in the 5 data sets of dif- 
ferent B (2.5 AiT, 5fiT, 7.5 AiT, lOfiT, and 12.5 ^fT) by 
fitting Eq. to the data. None of the fits showed a 
significant modulation. Limits on the amplitude were 
calculated performing a frequentist confidence level anal- 
ysis along the lines of [2l[ . The results of the fits and the 
corresponding limits are listed in Table |TT1 



B. Search for a Resonance 

In order to search for a resonance in the loss rate at 
the point B w B' , we averaged all normalised UCN 
counts for individual B field settings (thereby averag- 
ing out any remaining long term drifts) and plotted the 
results as a function of B (see Fig. [T]). A combined fit 
to the two data sets was performed using Eq. ([7]) with 
the following free parameters: two normalisation con- 
stants iV^g and N[^Q, the magnitude of B' , the angle 
P, and the oscillation time r„„/. The value for B' was 
constrained to lie in the region . . . 12.5 ^T as only in 
that region we would have unambiguous evidence for a 
possible resonance. The relevant, fitted parameters are 
B' = 11.4 ^T, l3 = 25.3 deg, and w = 21.9 s. The 



per degree of freedom (x^/dof = 17.86/17) is compara- 
ble to the one obtained by fitting a constant to the data 
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FIG. 2: (colour online) Contour plot of the minimal the 
point (£?', r„,j/). The solid line denotes the 95% C.L. contour 
line for an exclusion of t„„i . We evaluated a lower limit on 
T„„i at the minimum of this contour for B' between and 
12.5 fiT. 

(x^/dof = 22.72/21). There is therefore no evidence of 
a mirror magnetic field present at the site of the exper- 
iment and the data were used to set a limit on r„„' for 
mirror magnetic fields between and 12.5 f/T. To do so, 
the minimal x^ at the points {B' , r„„') was calculated 
by fitting the remaining free parameters ^^75, -^{501 ^^'^ 
P (see Fig. [5]). The 95% C.L. contour corresponds to 
X^ = 27.59, the 95% C.L. for a x^ distribution with 17 
degrees of freedom. Figure [2] also shows the loss of sen- 
sitivity to nn' oscillations for B' fields outside the range 
of applied magnetic fields. We evaluated a lower limit on 
the oscillation time as the minimal Tnn' on this contour 
for B' between and 12.5 fiT: 



Tnn' > 12.0 s (95% C.L. 



(11) 



The 0.1 /xT precision on individual non-zero B field val- 
ues leads in principle to a systematically improved limit. 
The improvement could not be quantified exactly, but it 
is estimated to be less than 1 s, and was not included in 
the result. Additionally, we improve our previous limit on 
Tnn' for negligible B' at the intercept of the exclusion con- 
tour line in Fig. [Hwith B' = 0: r„„/ > 141 s (95% C.L.). 
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